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Investigation of Surface Water Behavior During 

Glaze Ice Accretion' 

R. John Hansman Jr.* and Stephen R. Turnockf 
Massachusetts Institute of Technology , Cambridge , Massachusetts 

A series of experimental investigations that focused on isolating the primary factors that control the behavior of 
unfrozen surface water during glaze ice accretion were conducted. Detailed microvideo observations were made of 
glaze ice accretions on 2.54 cm diam cylinders in a closed-loop refrigerated wind tunnel. Distinct zones of surface 
water behavior were observed; a smooth wet zone in the stagnation region with a uniform water film, a rough zone 
where surface tension effects caused coalescence of surface water into stationary beads, and a zone where surface 
water ran back as rivulets. The location of the transition from the smooth to the rough zone was found to migrate 
towards the stagnation point with time. Comparative tests were conducted to study the effect of the substrate thermal 
and roughness properties on ice accretion. The importance of surface water behavior was evaluated by the addition 
of a surface tension reducing agent to the icing tunnel water supply, which significantly altered the accreted glaze 
ice shape. Measurements were made to determine the contact angle behavior of water droplets on ice. A simple 
multizone modification to current glaze ice accretion models was proposed to include the observed surface 
roughness behavior. 


Nomenclature 

A = surface area, m 2 

C, = specific heat capacity of ice, J/kg K 

C p = specific heat capacity of air, J/kg K 

C w = specific heat capacity of water, J/kg K 
D = diffusion coefficient of water vapor in air, m 2 /s 
F = force, kg m/s 2 

h = local convective heat-transfer coefficient, W/m 2 K 
k = thermal conductivity of air, W/m K 

L f — latent heat of fusion of water, J/kg 

L 3 = latent heat of sublimation of water, J/kg 

L v = latent heat of vaporization of water, J/kg 

M" = local mass flux/time, kg/m 2 s 

()" = local heat flux/time, W/m 2 

r = recovery factor, 0.875(-) 

T'surf = equilibrium surface temperature, °C 
= cloud temperature, °C 
/ = icing time, s 

Vao = freestream velocity, m/s 

W = cloud liquid water content, g/m 3 
P = local collection efficiency (-) 

AO = contact angle hysteresis, deg 
AT^ = cloud supercooling = T^{°C) 

0 = contact angle, deg 

0 a = advancing contact angle, deg 

0 r = receding contact angle, deg 

Pr.surf = saturated vapor density over surface, kg/m 3 

Pv , oo = saturated vapor density in cloud, kg/m 3 

a = surface tension, kg/m s 2 
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I. Introduction 

T HE behavior of unfrozen water on an accreting glaze ice 
surface can directly and indirectly influence the shape of 
the resulting ice accretion.' The surface transport behavior of 
the unfrozen water prior to freezing has a direct impact on the 
ice shape due to its effect of redistributing the impinging water 
mass. The surface water behavior also indirectly influences the 
ice accretion through its impact on surface roughness, which 
modifies the local convective heat transfer. The local heat 
transfer is the controlling factor in wet surface glaze ice accre- 
tion, where there is sufficient impinging water such that ice 
accretion rate is limited by the ability to remove latent heat of 
fusion from the surface. 

Current ice accretion models ignore the specific details of 
unfrozen surface water behavior during glaze accretion. The 
transport behavior and surface roughness are both treated in a 
simplistic or heuristic manner. This simplistic evaluation of 
surface roughness and transport behavior, which omits consid- 
eration of the surface physics, is considered to be a contribut- 
ing factor in the poor agreement between current glaze ice 
accretion models and expeimental results. 1 

This paper describes a series of experimental investigations 
focused on isolating the primary factors that control surface 
water behavior. The experimental investigations include mea- 
surement of the contact angle and resistance to motion of water 
on an ice surface, photographic observations of surface rough- 
ness during glaze ice accretion on cylinders, differential com- 
parisons of substrate thermal and surface properties on ice 
accretion, and the influence of surface tension on glaze ice 
accretion. The results indicate the importance of surface water 
behavior to the icing process in the glaze ice regime. A simple 
modification to existing glaze ice models is also proposed. 

II. Ice Accretion Mechanisms Influenced by Surface 
Water Behavior 

The primary ice accretion mechanisms that are influenced by 
the behavior of unfrozen surface water are the convective heat 
transfer, which is strongly related to the surface roughness, 2 
and the mass transfer at the surface, which is partly determined 
by the water runback characteristics of the surface. Both the 
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heat and the mass transfer will be discussed briefly in the con- 
text of ice accretion models to identify the importance of sur- 
face water behavior on the ice accretion process. 

Mass Transfer 

The principal modes of liquid water mass transfer on an 
accreting ice surface are shown schematically in Fig. 1, which 
depicts mass flux into and out of a control volume at the ice 
surface. The primary source of liquid water is the impingement 
of supercooled droplets. The impinging liquid water mass flux 
A/" mp is linearly related to the ambient liquid water content W, 
the local droplet collection efficiency 0, and the freestream 
velocity K*,: 

( 1 ) 

The double prime superscript is used to indicate that the 
quantity is defined per unit area of the icing surface. 

Liquid water may also enter the control volume through 
mass flow along the surface (M in in Fig. 1). Liquid water leaves 
the control volume through freezing at the ice water interface 
Ml rcezc , and surface flow out of the control volume A/ oul . Some 
mass also leaves the control volume due to evaporation, but 
this is generally small in the icing problem and will be neglected 
in the following analysis. The freezing mass flux Mf rccze is de- 
termined by the ability to convect the latent heat resulting from 
the water ice phase transition away from the surface. This will 
be discussed more fully in the following section. 

In steady state, the mass flux into and out of a control vol- 
ume with surface area A must balance 


gest two potential scenarios for surface water flow in the pres- 
ence of water beads. In one scenario, the water beads will grow 
until the increased heat transfer from the rough surface be- 
comes sufficient for the freezing mass flux to equal the imping- 
ing mass flux and any surface flow into the control volume. In 
this case, no water would leave the control volume, and Eq. (3) 
would reduce to 

A ■ A/f r ecze = A * A/'j' m p + A^ in (4) 

In the alternate scenario, the water beads would grow 
through coalescence and impingement until external dynamic 
pressure and shear forces caused the droplets to slide along or 
be torn from the surface. In this case, it is unlikely that the drop 
would move into the next control volume, and the assumption 
of surface flow continuity implicit in current mass balance 
models does not hold. This may explain experimental observa- 
tions of reduced total ice accretion at temperatures near 0°C. 5 

It should be noted that both scenarios described above may 
occur depending on the particular temperature and liquid wa- 
ter conditions. It should also be noted that uniform water films 
were also observed in the NASA experiments, particularly near 
the stagnation region. 

Heat Transfer 

The thermodynamic analysis presented in this paper for an 
accreting ice surface follows the earlier work by Olsen et al. 
and others, 6 7 and is commonly employed in current ice accre- 
tion models. 1 8 The principal modes of energy transfer associ- 
ated with the icing surface are depicted schematically in Fig. 3. 
Heat is added to the surface, primarily from the latent heat of 


A ■ A/; mp + M in = A ■ Mf recze + M out (2) 

For rime icing conditions at cold temperatures, the convec- 
tive heat transfer is sufficient to freeze all of the impinging mass 
flux. In this case, there is no surface flow, and Eq. (2) reduces 
to 

= A/ f reeze (3) 

In glaze icing conditions, there is insufficient heat transfer to 
freeze all of the impinging mass flux. In this case, the surface 
flow terms in Eq. (2) have to be considered. In current ice 
accretion models, all unfrozen water is assumed to run back 
out of the control volume, while unfrozen water from the next 
upstream element is assumed to flow into the control volume. 1 

The assumption of uniform runback is quite reasonable for 
icing conditions where a thin water film covers the glaze ice 
surface. However, recent NASA high-speed photographic 
studies indicate that other runback mechanisms may exist. 3 * 4 In 
some of the NASA studies, surface water was observed to 
coalesce into beads that remained stationary on the accreting 
ice surface. An example is shown in Fig. 2. The liquid beads 
were often surrounded by regions of otherwise dry ice surface, 
and no runback was observed. The NASA observations sug- 
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Fig. 2 Grazing angle photograph of a glaze ice surface (courtesy of 
William Olsen, NASA Lewis Research Center). 4 
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Fig. 3 Modes of energy transfer for an accreting ice surface. 
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fusion released as the droplets freeze, but also from aerody- 
namic heating and, to an even smaller extent, from the kinetic 
energy of the droplets impacting the surface. Heat is removed 
from the surface primarily by convection, and to a lesser degree 
by sublimation (when the surface is dry) or evaporation (when 
the surface is wet). In addition, heat is absorbed from the 
surface as the supercooled droplets impinge and warm to 0°C. 

The assumption of steady state requires that the rate at 
which energy is added to the control volume equals the rate at 
which it is removed, i.e., 

C^in “ Q. out ( 5) 

where £•„ and £" ut represent the energy added to and removed 
from the control volume per unit area per unit time. Equation 
(2) may be expanded into its component energy terms as 

^in ^ freezing aero heating droplet kinetic energy ( 6 ) 
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to the local convective heat-transfer coefficient. For wet sur- 
faces, the ice accretion rate is essentially limited only by the 
capability to transfer the latent heal of the freezing water away 
from the surface. 

The convective heat-transfer coefficient h has fcn;en empiri- 
cally observed to be strongly dependent on the ambient surface 
roughness. 2 Current ice accretion models either empirically 
generate local convective heat-transfer distributions to match 
natural ice shapes, 9 or assume a uniform roughness and use 
integral boundary-layer techniques to derive heat-transfer dis- 
tributions. 10 In these cases, the surface roughness is character- 
ized by an equivalent sand-grain roughness, which is adjusted 
to match experimental ice accretions. 

The surface water behavior is a primary factor in the evolu- 
tion of roughness on an accreting glaze ice surface. The surface 
roughness determines the heat-transfer behavior, which ulti- 
mately determines the ice accretion rate under glaze conditions. 
Understanding the generation of surface roughness by the mi- 
crophysical surface water behavior is, therefore, a key compo- 
nent to understanding the glaze ice accretion process. 


At steady state it is assumed that the ice surface achieves a 
locally uniform equilibrium temperature T surf . Conduction into 
the ice is assumed to be zero, and chordwise conduction be- 
tween adjacent control volumes is neglected. With these as- 
sumptions, the component heat terms of Eqs. (6) and (7) mav 
be written as 


@ freezing — ^freei c[Lf + C/(0°C — 7~ sur f)] 
6 : rh v - 
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Note that both of the primary heat dissipation terms ()" onv 
6evap contain the local convective heat-transfer coefficient 
h. Also note that the term that dominates the heat flux 

into the control volume contains the freezing mass flux Aff rc€ „. 
This implies that the local ice accretion rate is strongly coupled 


Fig. 4 Example of contact angle 
hysteresis in a droplet on a vertical 
wall. 



III. Microphysics of Surface Water Behavior 

Forces Influencing Surface Water Behavior 

The microphysical behavior of water on an aircraft surface is 
controlled primarily by the relative strength of the surface ten- 
sion, and by aerodynamic and body forces. Surface tension 
forces tend to minimize the surface area of the fluid, causing 
the water to coalesce into beads or rivulets. They also act to 
oppose motion of fluid along the surface. Body forces such as 
gravity or centripetal acceleration act on the entire fluid bulk. 
Aerodynamic forces result from pressure gradients and shear 
stress at the water-air interface. 

Contact Angle 

An important parameter in surface fluid behavior is the con- 
tact angle 9. It is defined as the angle the fluid-vapor interface 
makes with an underlying substrate. The nominal contact angle 
is a property of the specific gas-liquid-solid combination, and 
may vary with temperature. The wettability, or tendency of a 
fluid to spread on a particular surface, is inversely related to the 
contact angle d. In general, for contact angle values of less 
than 10 deg, the surface is considered wettable. 

If the fluid is subject to an external force such as gravity, the 
contact angle may vary, depending on the direction of the force 
relative to the contact line. This hysteresis effect can be visual- 
ized in the schematic view of a drop sliding down a vertical 
surface, shown in Fig. 4. On the lower surface of the drop the 
advancing line of contact causes an increase in the contact 
angle to its maximum allowable value 0 a . Conversely, the re- 
ceding edge of the drop remains at the lowest allowable value 
of the contact angle 6 r . The contact angle hysteresis A6 is the 
difference between the advancing and receding contact angles: 

AO = 6 a - 6 r (9) 

The contact angle hysteresis AO is a property of the specific 
gas-hquid-solid combination, and will tend to increase with 
surface roughness. The resistance to motion of the liquid-solid 
line of contact can be related to the value of the nominal con- 
tact angle and its hysteresis. 11 In general, increasing values of 0 
and AO imply increased resistance to motion. 

Modes of Surface Water Behavior 

Water on an aircraft surface tends to behave in one of three 
distinct modes depending on the water flow rate, aerodynamic 
forces, and the wettability of the surface. For high flow rates 
and wettable surfaces, the water will tend to coat the entire 
surface with a uniform film. For lower flow rates, or higher 
contact angles, there is insufficient water to coat the entire 
surface. Once drying begins, surface tension forces will tend to 
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coalesce the available water into either rivulets or water beads. 

Water beads are generally associated with low flow rates or 
unwettable surfaces such as wax. Equation ( 10) represents the 
overall force balance for a small stationary water bead. As long 
as the contact line force is greater than the sum of the aerody- 
namic and body forces, the bead will remain stationary. 

^contact ^ ^aero "1” ^"body ( 

The contact line force F c oma ci IS the tota ^ force resulting from 
the contact line resistance to motion described in the preceding 
section. Both the contact line force and the aerodynamic force 
are roughly proportional to the beads surface area, whereas the 
body force is proportional to the bead volume. This implies 
that, as the bead size is increased, there will be maximum 
stationary bead volume for a given set of surface, flow, and 
gravitational conditions above which the bead will move. 

IV. Measurement of Contact Angle Between Ice and 
Water 

As discussed in Sec. Ill, the contact angle 9 and hysteresis A 9 
are important parameters in the behavior of water on an ice 
surface. There is, however, very little data available on 6 and 
A 9 for water on ice. This is partly due to experimental 
difficulties. Most investigations of the surface properties of su- 
percooled water have concentrated on measurements of surface 
tension a. 12 In order to obtain preliminary working values of 9 
and Ad for water on ice, a series of simple experiments were 
conducted. 

Experimental Setup 

The apparatus employed in these investigations is shown 
schematically in Fig. 5. A smooth layer of ice, approximately 
8 mm thick, was formed from distilled water on a metal plate 
that could be oriented at a variety of angles with respect to the 
horizontal. Precooled droplets of distilled water were placed on 
the ice surface by a syringe, and their shapes were recorded by 
a CCD microvideo camera oriented to view the ice surface at a 
grazing angle. By use of a dark background and a diffuse light 
source, a sharp high-contrast image of the droplet at the ice- 
water interface could be obtained. This enabled contact angles 
to be measured to an accuracy of approximately ± 5 deg. The 
ice surface temperature was monitored by an iron-constantan 
thermocouple mounted in direct contact with the ice surface. 
The output of the thermocouple was presented on a digital 
display within the field of view of the video camera. This pro- 
duced a simultaneous record of contact angle and surface tem- 
perature valid to approximately ± 1°C. 

Measurement of the contact angle at various temperatures 
was obtained by initially cooling the iced rrxetal plate in a cold 
box, to a temperature below the target measurement tempera- 
ture. The plate was then removed from the cold box and 
mounted horizontally within the field of view of the camera. As 
the plate warmed to the target temperature, droplets were 
placed on the surface and their shapes were recorded. Although 
the droplets would begin to freeze shortly after placement, the 
initial contact angle was preserved in the frozen drop, and any 
transient effects could be observed in the video record. In this 
manner, values of the contact angle of water on ice could be 
obtained for subfreezing temperatures. 

The contact angle hysteresis measurements were made utiliz- 
ing a similar procedure. However, the iced plate was set to an 
angle of 30 deg with respect to the horizontal. Each droplet’s 
volume was increased by syringe injection until motion began; 
at this instant the advancing and receding angles and the sur- 
face temperature were measured from the video recording. 
Results 

The observed dependence of contact, angle 9 and con tael 
angle hysteresis AO with temperature are shown in Figs. 6 and 
7, respectively. Because of the low magnitude, it was not possi- 
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Fig. 5 Schematic of contact angle measurement setup. 



ble to obtain accurate contact angle measurements at tempera- 
tures above — 4°C. 

However, at the freezing point (0°C), ice and water are in 
equilibrium so that the ice surface must be perfectly wettable, 
implying that both 0 and AO must be zero. These points are 
included in the contact angle and contact angle hysteresis plots. 
They are consistent with the experimentally observed values 
and allow interpolation between — 4°C and 0°C. 

Both the contact angle and hysteresis exhibit a strong non- 
linear variation with temperature, particularly in the vicinity of 
the freezing point (0°C). These results imply that the wettabil- 
ity of the ice surface decreases greatly as the ice surface cools 
below freezing. For warm ice surfaces (near 0°C) water will 
tend to spread, and the ice-water contact line will have a low 
resistance to motion. For colder ice temperatures, water will 
tend to bead into droplets, which require a higher force to 
initiate fluid motion. 

The strong contact angle temperature dependance indicates 
the potential importance of thermal gradients to the develop- 
ment of ice surface roughness. Small variations in surface tem- 
perature can restrict the mobility of water. This effect is 
thought to be the cause of the stable nature of the surface water 
beads observed by Olsen and Walker 3,4 and shown in Fig. 2. 
The dry ice surface around the beads can be cooled to temper- 
atures well below 0°C by convective heat transfer. This cold dry 
surface, therefore, imposes a barrier to water flow away from 
the bead. Any impinging water that strikes the dry surface will 
quickly be frozen due to the low temperature. Any cloud 
droplets that strike the water will be trapped within the bead, 
while glaciation occurs at the lower surface of the bead. In this 
manner, the ice thickness can increase, while the bead is ob- 
served to remain stationary. 

V. Experimental Observations of Ice Surface Roughness 
on Cylinders 

Experiment Setup 

Detailed photographic observations of the behavior of sur- 
face water and formation of ice roughness on cylinders during 
glaze ice accretion were made in the Data Products of New 
England Icing Wind Tunnel. The experimental setup is shown 
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in Fig. 8. The tunnel was a closed-loop refrigerated system with 
a 6-in. -square plexiglass test section. The test article was a 
cylinder that horizontally spanned the test section. Scale refer- 
ence for the photographic studies was provided by a grid 
mounted on a thin splitter plate at the midplane of the test 
section. An 8 mm CCD microvideo (camera A) with a macro- 
lens for magnification was used to obtain a grazing angle view 
of the ice accretion. The camera was focused at the stagnation 
region near the center of the test section. An additional video 
camera (camera B) was mounted above and slightly upstream 
of the cylinder to provide a view of the ice accretion looking 
normal to the cylinder surface. Secondary lighting was pro- 
vided to obtain a clear video record of the accreting ice surface. 

The tunnel and icing cloud parameters were set to yield a 
variety of glaze ice conditions. The air velocity was 150 knots 
for all tests. A range of freestream air temperatures between 
—9 and — 4°C was used. The centerline liquid water content W 
was varied from 0.7 to 1.2 g/m 3 . The icing cloud for the condi- 
tions had a nominal mean volumetric diameter (MVD) of 30 p. 
The standard exposure time was 3 min. The tunnel had previ- 
ously been calibrated using an indirect method. There is, there- 
fore, some uncertainty in the liquid water content and MVD 
values. 

The liquid water content was observed to vary somewhat 
across the test section with a reasonably constant valve in the 
central 2 in. as measured by rime ice accretion calibration runs. 
Photographic observations were, therefore, focused on the cen- 
terline region where the ice deposit was uniform and the liquid 
water content calibration was valid. 




Results 

Figure 9 shows ai representative final glaze ice shape for an » 
exposure temperature ol — 4.5 C and a liquid water contact of 
1.0 g/m 1 . During the accretion, the following three distinct 
types of ice surface behavior were observed, each having a 
characteristic roughness and identifiable boundaries. 

1) Smooth Zone . Close to the stagnation point during expo- 
sure to the icing cloud, the surface was observed by light reflec- 
tion techniques to be uniformly wet with a thin film of water at 
warm temperatures. The surface in this regime was smooth, 
with no distinct roughness. The ice was translucent within this 
zone. 

2) Rough Zone. At some point downstream, there was a 
sudden transition to a significantly rougher surface. Within this 
zone, there appeared to be insufficient water to maintain a 
uniform film. Surface tension forces dominated the water sur- 
face behavior. Runback did not occur; rather, the water tended 
to coalesce into the water beads first observed by Olsen and 
Walker. 3 - 4 The scale length of the roughness was typically on 
the order of 1 mm. The transition between the rough and the 
smooth zones can be clearly seen in the grazing angle 
photograph shown in Fig. 2. 

Inasmuch as there was a distinct boundary between the 
smooth and rough zones, this position could be easily identified 
on the grazing angle video recording. The angular position of 
this boundary is plotted in Fig. 10 for cold ( -9°C) and warm 
( -4.5' C) conditions at a liquid water content of 1.0 g/m 3 . The 
boundary started at approximately 50 deg and propagated 
rapidly towards the stagnation region. The repeatable nature 
of the smooth-rough transition’s propagation toward the stag- 
nation region implies a clear underlying physical mechanism 



Fig. 9 Typical final glaze ice shape showing distinct roughness zones, 
~4.5°C, 1 .0 g/m \ 



Fig. 1 0 Angular position vs time of smooth-rough transition location for 
two air temperatures. 
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for the transition between the surface water behavior in the 
smooth and rough zones. 

The ice accretion rate was observed to be enhanced in the 
rough zone as compared with the smooth zone. This may be 
seen in the ice profiles of Fig. 1 1 taken from the video at 30, 90, 
and 150 s after exposure to the icing cloud. The enhanced ac- 
cretion rate is thought to be due to increased heat transfer 
resulting from the greater surface roughness in this zone. In all 
of the cases observed in this study, the ice horns characteristic 
of glaze ice accretion were found within the rough zone. ^ 

3) Runback Zone. At warm temperatures a third zone was 
observed aft of the rough zone. This region was characterized 
by areas of ice interspersed with uniced surface. This ice was 
observed to form during an initial transition period after cloud 
exposure. The ice was translucent, and quite often frozen 
rivulets could be discerned. In warm conditions and high liquid 
water contents, the surface water was observed to initially run- 
back and then stagnate at the point of flow separation. This 
water then slowly froze as rivulets or as large coalesced water 
cells. Once ice began to form in the upstream rough zone, no 
additional surface water was supplied to the runback zone, and 
the ice surface remained constant. 

VI. Observations of the Effect of Substrate Properties 
on Surface Water Behavior 
Experimental Setup 

Two properties of the initial uniced surface that can affect 
the surface water behavior are 1) thermal characteristics (i.e., 
conductivity, specific heat capacity), and 2) roughness of the 
uniced surface. By exposing two test articles (identical except 
for either their thermal characteristics or initial surface rough- 
ness) to the same icing cloud, it was possible to assess the 
individual importance of each parameter. Comparative tests of 
this type were carried out in the icing wind-tunnel setup de- 
scribed in Sec. V. For these experiments, the test article was 
composed of two 3-in.-long, l-in.-diam cylinders of different 
material composition, separated by a splitter plate. As shown 
in Fig. 8, each half of the test article was viewed by a grazing 
angle microvideo camera (cameras A and C). A third camera 
(B) was mounted above and upstream of the article to simulta- 
neously record the differences between the two cylinders. 

Two cylinders were used in the thermal comparison experi- 
ments! a copper tube with a 1/16-in. wall thickness to investi- 
gate fast thermal response, and a solid plexiglass rod to 
investigate slow thermal response. To remove the possible in- 
fluence of surface chemistry effects, each cylinder was covered 
with a single coat of acrylic paint. 

For the roughness experiment, two cylinders were used. Both 
were manufactured from solid aluminum rods. One was ex- 
tremely smooth with a polished finish obtained using ‘0000’ 
emery paper. The rough cylinder had a repeatable surface pat- 
tern that was produced by knurling the cylinder on a lathe. The 
knurling process produces a pattern of trapezoidal surface ele- 
ments. For the cylinder used in these experiments, the elements 
had a height of 0.8 mm, a width of l mm, and a length of 

2.7 mm. 

Results 

Thermal Comparison 

Figure 12 shows a typical comparison of the final ice shapes 
obtained after 3 min exposure for the copper tube and the 
plexiglass rod. This comparison was run at a temperature of 
-5.5"C and a liquid water content of 0.95 g/m 3 . Although 
there is not a great deal of difference between the final accre- 
tions, the glaze ice horns on the copper cylinder were more 
sharply defined than on the plexiglass, which had a relatively 
flat front surface and was slightly thicker in the stagnation 
region. The angular position of the rough-smooth transition 
boundary is shown as a function of exposure time for these 



Fig. 1 1 Time development of glaze ice shapes, T = — 4.5°C, 

W = 1.0 g/m 3 . 


y Plexiglass 



Fig. 12 Typical final ice shape comparison for copper and plexiglass 
cylinders, T « -5.5°C, W = 0.95 g/m 3 . 



Time (secs) 

Fig. 13 Angular position of smooth-rough transition vs time for copper 
and plexiglass cylinders. 


accretions in Fig. 13. For the plexiglass cylinder, the smooth 
zone shrinks much more rapidly than for the copper. The ice 
surface on the plexiglass cylinder became uniformly rough after 
45 s of exposure. The enhanced roughness in the stagnation 
region, after this time, accounts for the thicker ice in the stag- 
nation region and the relatively flat front ice surface observed 
in Fig. 12. 
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The water surface behavior after initial exposure to the icing 
cloud was observed to differ between the copper and plexiglass 
test articles. For the copper cylinder, there was relatively little 
surface water flow prior toTreezing, and the initial icc surface 
was relatively smooth. For the plexiglass, droplet coalescence 
and runback were observed to occur prior to the initial ice 
formation. This resulted in a significantly rougher initial ice 
surface, which is thought to be the cause of the difference in 
accretion behavior observed between the plexiglass and copper 
test articles. 

In general, the differences in ice accretion behavior due to 
substrate thermal properties were observed to be most signifi- 
cant during the initial phase of ice accretion. As the accretion 
grows, the effect of initial conditions begins to wash out. After 
extended icing encounters, the accretion will tend toward a 
shape that is controlled by the environmental parameters and 
is independent of substrate properties. 

Vniced Surface Roughness 

A typical comparison of final ice profiles for the polished 
(smooth) and knurled (rough) cylinders is shown in Fig. 14. 
The ice horns on the knurled cylinder were broader and further 
off the stagnation line than the horns on the polished cylinder. 
This is thought to be due to enhanced heat transfer in the horn 
regions during the initial phase of the ice accretion. A slight but 
repeatable asymmetry was observed in the polished ice growth. 
This is thought to be due to the smooth surface on the polished 
cylinder, which allowed gravitational forces to influence the 
initial surface mass flux, resulting in a larger ice horn on the 
lower surface. 

The zones of smooth and rough ice growth discussed in Sec. 
V were observed in both the polished and knurled cases. Figure 
15 plots the angular position of the rough-smooth transition vs 
time. There is some difference between the two cases, with the 
smooth zone being smaller for the knurled cylinder after ap- 
proximately 80 s exposure time. An important observation was 
that the trapezoidal element pattern of the knurled cylinder 
was clearly apparent in the ice accretion within the rough zone. 
This implies that the roughness geometry propagates in the ice 
accretion within the rough zone. In the absence of uniced sur- 
face roughness, water beading mechanisms appear to be the 
source elements for the roughness geometry. 

VII. Surface Tension Effects on Ice Accretion 
Experimental Setup 

The importance of surface water behavior was further inves- 
tigated by varying the surface tension of the water and studying 
the effect on the resulting ice accretion. The water surface ten- 
sion was varied by the addition of a surfactant (Photoflo 200), 
which reduced the surface tension by approximately a factor of 
2, while leaving its bulk properties (density, freezing tempera- 
ture) unchanged. 

The icing wind-tunnel setup was similar to that discussed in 
Sec. V. Photoflo was added at a dilution of 1:200 to the icing 
cloud water supply. Icing tunnel conditions were matched as 
closely as possible between runs carried out with and without 
the surfactant added. Because the spray system was not specifi- 
cally calibrated with Photoflo, some uncertainty exists on liq- 
uid water content and MVD values. Based on splitter plate ice 
accretion, this uncertainty is estimated to be on the order of 
5%. Photoflo runs were conducted at a liquid water content of 
1 g/m and for temperatures between -9 and -4°C 
Results 

Final glaze ice profile comparisons between a normal icing 
cloud and the Photoflo modified cloud exposures are shown in 
Fig. 16 for the copper cylinder at -4.5’C. The addition of 
Photoflo resulted in an opaque white icc with a smaller scale 
surface roughness and a significantly different profile. The 
horns were more pronounced and closer to the stagnation 
point. At the lower temperatures, the glaze ice horn has bc- 



Fig. 14 Typical final ice shape comparison for polished and knurled 
cylinders, T » -4.5*C, W = 1.0 g/m 3 . 




Fig. 16 Typical final ice shape comparison for cloud 
and without surfactant, T = -4.5°C, W = 0.95 g/m\ 


water supply with 
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come so pronounced that it no longer collected all of the im- 
pinging water. A second horn was, therefore, able to develop 
behind it. 

The different and unusual ice accretions obtained at reduced 
values of surface tension indicate that surface water behavior 
is, indeed, an important factor in the glaze ice accretion process 
and must be considered in physically realistic models. 

VIII. Implications for Ice Accretion Modeling 

Based on the observations of surface water behavior during ^ 
glaze ice accretion, a relatively simple modification to the exist- 
ing ice accretion model is proposed, which may improve the 
current model's accuracy within the glaze ice regime. In this 
proposed “multi-zone" model, the accreting ice surface is di- 
vided into two or more discrete zones that have varying surface 
water behavior and surface roughness. This is in contrast to 
current techniques, which assume that the surface has uniform 
roughness and surface water runback. 

In the simplest version of the model, the surface is divided 
into two zones. There is a smooth wet zone centered about the 
stagnation region where thin film runback occurs, and the ex- 
isting Messinger-type models 6 appear to be valid. In this region 
the heat transfer would be that for a smooth surface. The 
remaining ice surface would consist of a “rough" zone where 
surface tension effects dominate the surface water behavior, 
and the characteristic water beads appear. The heat transfer 
within this region would be adjusted by incorporation 
of roughness correction factors appropriate to the bead 
dimensions. 

Two primary areas need to be addressed prior to implemen- 
tation of a two-zone model. They are: 1) determining the 
smooth-rough transition locations, and 2) determining the ap- 
propriate value of equivalent sand-grain roughness to be used 
in calculating the rough zone heat-transfer coefficient. 

The work described in Sec. V indicates that there are distinct 
physical mechanisms that control the position of the smooth- 
rough transition. However, further work is necessary to para- 
metrically define the location of the smooth-rough transition. 

In addition, the factors that control the generation of rough- 
ness elements within the rough zone need to be investigated in 
the context of their impact on the local convective heat-transfer 
coefficient. 

IX. Conclusions 

The investigation of surface water behavior on an accreting 
glaze ice surface has shown the following: 

1) Surfaces have been observed to have distinct zones of 
surface water behavior. They include a smooth wet zone in the 
stagnation region where uniform film runback occurs, a rough 
zone where surface tension causes coalescence of the surface 
water into beads and no runback occurs, and a runback zone 
where surface water runs back as rivulets. 

2) The location of the transition point between the smooth 
and rough zone was observed to propagate with time toward 
the stagnation point. 

3) The freezing of the coalesced water beads in the rough 
zone generates a characteristic rough surface that enhances 
heat transfer. 


4) Initial surface roughness patterns were observed to prop- 
agate in the accreted ice within the “rough zone." 

5) Large variations in ice accretions were observed when the 
surface tension of the water was changed. This illustrates the 
importance of surface water behavior. 

6) Measurement of contact angle and contact angle hys- 
teresis for water on an ice surface showed a strong variation 
with the ice surface temperature. This temperature dependence 
indicates that thermal gradients along the ice surface may be 
important to surface water behavior and to the generation of 
roughness on an accreting ice surface. 

7) Based on this investigation, a simple multiple-zone mod- 
ification to the current glaze ice accretion model is proposed. 
The model incorporates discrete zones of surface water behav- 
ior each with a characteristic surface roughness. 
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